Nuclear symmetry energy and the r-mode instability of neutron stars 
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We analyze the role of the symmetry energy slope parameter L on the r-mode instability of neutron stars. Our 
study is performed using both microscopic and phenomenological approaches of the nuclear equation of state. 
The microscopic ones include the Brueckner-Hartree-Fock approximation, the well known variational equation 
of state of Akmal, Pandharipande and Ravenhall, and a parametrization of recent Auxiliary Field Diffusion 
Monte Carlo calculations. For the phenomenological approaches, we use several Skyrme forces and relativisic 
mean field models. Our results show that the r-mode instability region is smaller for those models which give 
larger values of L. The reason is that both bulk {£,) and shear (rj) viscosities increase with L and, therefore, the 
damping of the mode is more efficient for the models with larger L. We show also that the dependence of both 
viscosities on L can be described at each density by simple power-laws of the type £, = A^lPi and 7] = Aj^L^'' . 
Using the measured spin frequency and the estimated core temperature of the pulsar in the low-mass X-ray 
binary 4U 1608-52, we conclude that observational data seem to favor values of L larger than --^ 50 MeV if this 
object is assumed to be outside the instability region, its radius is in the range 11.5 — 12(11.5 — 13) km, and its 
mass 1.4Mq)(2Mqj). Outside this range it is not possible to draw any conclusion on L from this pulsar. 
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I. INTRODUCTION 



It is well known that the absolute upper limit on the rota- 
tional frequency of a neutron star is set by its Kepler frequency 
^hinjier, above which matter is ejected from the star's equator 
|lll|211- However, a neutron star may be unstable against some 
perturbations which prevent it from reaching rotational fre- 
quencies as high as H^epien setting, therefore, a more strin- 
gent limit on its rotation |3]. Many different instabilities 
can operate in a neutron star Among them, the so-called r- 
mode instability, a toroidal mode of oscillation whose restor- 
ing force is the Coriolis one, is particularly interesting and, 
since its discovery by Andersson, Friedman and Morsink a 
few years ago |4, "s*], its study has received a lot of atten- 
tion |6-19]. This oscillation mode leads to the emission of 
gravitational waves in hot and rapidly rotating neutron stars 
through the Chandrasekhar-Friedman-Schutz (CFS) mecha- 
nism ll20l l2l[] . and it is generally unstable at all rotational 
frequencies [4]. Gravitational radiation makes an r-mode to 
grow, whereas viscosity stabilizes it. Therefore, an r-mode is 
unstable if the gravitational radiation driving time is shorter 
than the damping time scale due to viscous processes. In this 
case, a rapidly rotating neutron star could transfer a signifi- 
cant fraction of its rotational energy and angular momentum to 
the emitted gravitational waves. These waves, potentially de- 
tectable, could provide invaluable information on the internal 
structure of the star, and constraints on the nuclear equation of 
state (EoS) |22]. In particular, r-modes can help to restrict the 
density dependence of the symmetry energy, £■„.„, (p), a cru- 
cial ingredient of the nuclear EoS needed to understant many 
important properties of isospin-rich nuclei and neutron stars 
II23I425I1 . The value of the symmetry energy at saturation den- 
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sity po is more or less well established (po) ^ 30 MeV), 
and its behaviour below po is now much better known ll26ll . 
However, above po, Esym{p) is not well constrained yet, and 
the predictions from different approaches strongly diverge. 
The r-modes can provide information on Esym (p ) complemen- 
tary to the one obtained from: (i) the analysis of data of giant 
1 27] and pygmy f2§] resonances, (ii) isobaric analog states 
|29], (iii) isospin diffusion measurements |30], (iv) isoscaling 
1 31], (v) meson production in heavy ion collisions ll32i l33ll . 
(vi) measurements of the neutron skin thickness in heavy nu- 
clei ]34-38], (vii) the characterization of the core-crust tran- 
sition in neutron stars lf39l - l42ll . (viii) the analysis of power- 
law correlations, such as the relation between the radius of a 
neutron star and the EoS ll43ll . or (ix) the novel constraints re- 
cently reported by Steiner and Gandolfi ]57] on the basis of 
neutron star mass and radius measurements, driven partially 
by the strong correlation between the symmetry energy and 
its derivative obtained in Quantum Monte Carlo calculations 
of neutron matter. 

In this work we want to study the role of the symmetry en- 
ergy slope parameter L — 3 po[dEsym{p)/dp]po on the r-mode 
instability. A similar study has been recently done by Wen, 
Newton and Li in Ref. ]45] using a simple model for the EoS 
that consistently describes the crust-core transition density. 
Assuming that the main dissipation mechanism of the r-modes 
is due to the electron-electron scattering at the crust-core 
boundary, and using the estimated core temperature of several 
low-mass X-ray binaries (LMXB), these authors conclude that 
neutron stars are stabilized against r-mode oscillations if L is 
smaller than ^ 65 MeV. In our work, we employ the micro- 
scopic Brueckner-Hartree-Fock (BHF) approach and several 
phenomenological Skyrme forces and relativistic mean field 
models to describe the neutron star matter EoS. In addition, 
we also use two other microscopically based EoS: the well 
known variational Akmal-Pandharipande-Ravenhall (APR) 
EoS [i461 . and the parametrization of the recent Auxiliry Field 
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Diffusion Monte Carlo (AFDMC) calculation of Gandolfi et 
al, given in Ref. ll47ll . We consider both the bulk (^) and 
shear (t]) as the main dissipative mechanisms of r-modes, in- 
cluding in the calculation of <^ the contribution of the modified 
and direct electron and muon Urea processes and, in that of r\, 
the contribution of neutron and electron scattering. 

The paper is organized in the following way. In Sec.HIlwe 
briefly review the BHF approach of nuclear matter and pro- 
vide some details on the Skyrme forces and the relativistic 
models considered. Few words about the APR and AFDMC 
EoS are also said. Section Hill is devoted to the calculation of 
bulk and shear viscosities. The dissipative time scales of the 
r-modes are presented in Sec. |IV| whereas the r-mode insta- 
bility region is determined in Sec. [V] Finally, a summary of 
our results is given in Sec. |VT1 



n. NUCLEAR EQUATION OF STATE MODELS 

The BHF approach is the lowest order of the Brueckner- 
Bethe-Goldston (BBG) many -body theory In this the- 
ory, the ground state energy of nuclear matter is evaluated 
in terms of the so-called hole-line expansion, where the per- 
turbative diagrams are grouped according to the number of 
independent hole-lines. The expansion is derived by means 
of the in-medium two-body scattering G matrix. The G ma- 
trix, which takes into account the effect of the Pauli princi- 
ple on the scattered particles and the in-medium potential felt 
by each nucleon, has a regular behavior even for short-range 
repulsions, and it describes the effective interaction between 
two nucleons in the presence of a surrounding medium. In 
the BHF approach, the energy is given by the sum of only 
two-hole-line diagrams including the effect of two-body cor- 
relations through the G matrix. It has been shown by Song 
et al. [49] that the contribution to the energy from three- 
hole-line diagrams (which account for the effect of three-body 
correlations) is minimized when the so-called continuous pre- 
scription 1 50] is adopted for the in-medium potential, which 
is a strong indication of the convergence of the hole-line ex- 
pansion. We adopt this prescription in our BHF calculation 
which is done using the Argonne VI 8 nucleon-nucleon poten- 
tial ifsill supplemented with a three-body force of Urbana type 
ll52ll . which for the use in the BHF calculation was reduced 
to a two-body density dependent force by averaging over the 
spatial, spin, and isospin coordinates of the third nucleon in 
the medium [53]. This three-body force contains two parame- 
ters that are fixed by requiring that the BHF calculation repro- 
duces the energy and saturation density of symmetric nuclear 
matter (see Refs. ll54l - l56il for a recent analysis of the use of 
three-body forces in nuclear and neutron matter). The inter- 
ested reader is referred to Ref. ll48ll for an extensive review of 
the BBG many-body theory, and to Ref. fST] for the specific 
details of our BHF calculation of isospin asymmetric nuclear 
matter Regarding the other two microscopic approaches used 
in this work, we note here that the APR and the AFDMC EoS 
have been implemented using, respectively, the parametriza- 
tions given by Heiselberg and Hjorth-Jensen in Ref. 15811 . and 
by Gandolfi et al., in Ref. li47il . 



Phenomenological approaches, either relativistic or non- 
relativistic, are based on effective interactions that are fre- 
quently built to reproduce the properties of nuclei. Skyrme 
interactions |59] and relativistic mean field models |60] are 
among the most commonly used ones. Many of such interac- 
tions are built to describe systems close to the isospin symmet- 
ric case, therefore, predictions at higher isospin asymmetries 
should be taken with care. Most of the Skyrme forces are, 
by construction, well behaved close to the saturation density 
and moderate isospin asymmetries. Nevertheless, only certain 
combinations of the parameters of these forces are well de- 
termined empirically. Consequently, there is a proliferation of 
different Skrme interactions that produce a similar equation of 
state for symmetric nuclear matter, but predict a very differ- 
ent one for pure neutron matter Few years ago. Stone et al., 
ll61i] tested extensively and systematically the capabilities of 
almost 90 existing Skyrme parametrizations to provide good 
neutron-star candidates. They found that only 27 of these 
parametrizations passed the restrictive tests they imposed, the 
key property being the density dependence of the symmetry 
energy. These forces are SLyO-SLylO \62i and SLy230a f&M 
of the Lyon group, SkIl-SkI5 |64] and SkI6 |65] of the Ski 
family, Rs and Gs [|66j], SGI ||67t], SkMP iH, SkO and SkO' 
[69], SkT4 and SkT5 fj^, and the early SV [TT]. The re- 
sults for the Skyrme forces shown in this work are restricted 
to these 27 parametrizations. We should mention, however, 
that more stringent constraints to the Skyrme forces have been 
very recently presented by Dutra et al, in Ref. [72]. These au- 
thors have examined the suitability of 240 Skyrme interactions 
with respect to eleven macroscopic constraints derived mainly 
from experimental data and the empirical properties of sym- 
metric nuclear matter at and close to saturation. They have 
found that only 5 of the 240 forces analyzed satisfy all the 
constraints imposed. We note that none of the 27 parametriza- 
tions used in this work corresponds to any of these 5 forces. 

Relativistic mean field models are based on effective La- 
grangian densities where the nucleon-nucleon interaction is 
described in terms of meson exchanges. In this work 
we consider two types of them: models with constant 
meson-nucleon couplings described by the Lagrangian den- 
sity of the nonlinear Walecka models (NLWM), and models 
with density-dependent couplings (hereafter called density- 
dependent hadronic models (DDHM)). In particular, within 
the first type, we consider the models GMl and GM3 JtU, 
TMl Q, NL3 and NL3-II and NL-SH For 
the DDHM models, we consider the models DDMEl and 
DDME2 ill, TW99 [78], and the models PKl, PKIR and 
PKDD of the Peking group ^791] . 

We finish this section by mentioning a few limitations of the 
models considered. Firstly, hyperons or other exotic degrees 
of freedom have not been considered in this work, although 
they are expected to appear in the inner core of neutron stars at 
few times saturation density. The reader should also note that 
causality is not always guaranteed for the BHF and AFDMC 
approaches, and the Skyrme forces. That is not surprising, 
since these models are non-relativistic. This is not the case, 
however, of the Heiselberg and Hjorth-Jensen parametriza- 
tion L58i] of the APR EoS where the EoS is softened at higher 
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FIG. 1: (Color online) Bulk (panel a) and shear (panel b) viscosities for the BHF calculation of nonsuperfluid j3-stable npefl matter as a 
function of the density for T = 10^ K and (o = 10* s^'. Contributions to the bulk viscosity from MURCA and DURCA processes involving 
electrons and muons are considered, as well as the contributions to the shear viscosity from neutron and electron scattering. 



densities to obey causality. Finally, the nuclear EoS has been 
obtained at zero temperature for each model. Temperature 
enters in our calculation only through the bulk and shear vis- 
cosities presented in the following section. 



III. BULK AND SHEAR VISCOSITIES 

Bulk and shear viscosities are usually considered the main 
dissipation mechanisms of r- and other pulsation modes in 
neutron stars. Bulk viscosity is the dominant one at high tem- 
peratures (T > 10^ K) and, therefore, it is important for hot 
young neutron stars. It is produced when the pulsation modes 
induce variations in the pressure and the density that drive the 
star away from j3 -equilibrium. As a result, energy is dissi- 
pated as the weak interaction tries to re-establish the equilib- 
rium. In the absence of hyperons or other exotic components, 
the bulk viscosity of neutron star matter is mainly determined 
by the reactions of modified Urea (MURCA) 

N + n^N + p + l + Vi , N + p + l^N + n + Vi 

{N^n,p , l^e-,11-) , (1) 



and direct Urea (DURCA) processes 

n ^ p + 1 + Vi , p + l ^n + Vi 



(2) 



the second one allowed only when the proton fraction Xp ex- 
ceeds a critical value xdurca Modified and 
direct Urea contributions to the bulk viscosity of nonsuper- 
fluid and superfluid jS -stabl e np e and npejJ. matter have been 
studied by several authors [8 1 In this work we assume 
that the neutron star interior is made only of neutrons, pro- 
tons, electrons and muons in a normal fluid state, and follow 
the work of Haensel et al. jsillssll to evaluate the bulk viscos- 
ity. According to these authors, since the Urea reaction rates 
are much smaller than the typical values of the frequency of 
the r-modes (in the range (O ^ 10^ — 10"* s^'), the total bulk 



viscosity can be simply written as a sum of the partial bulk 
viscosities associated with each (modified and direct) Urea 
process. 
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In the above expression co is the frequency of the pulsation 
mode, P is the pressure, ~ n„ + rip is the total baryon 
number density, Xi = rii/nt, is the electron or muon frac- 
tion, rji = pL„ — pLp — pLi, with pLi the chemical potential of 
the species /, and A^;/ and A/ determine the difference of the 
rates of the direct and inverse reactions of a given Urea pro- 
cess: — Tni — ^^Nili for the MURCA processes, and 
T/ - T/ = -A/Tj, for the DURCA ones. Note that when the 
system is in chemical equilibrium rji = 0. Note also that 
the quantities dP/dX; and drji/drii, depend on the particu- 
lar choice for the equation of state. The interested reader is 
referred to the original work of Haensel et al. llsiilssll for de- 
tails on the derivation of the specific expressions for the bulk 
viscosity employed here (see in particular Eq. (35)ofRef. llsill 
and Eqs. (18-21) of Ref. iH). 

Shear viscosity rj is the main viscous dissipation at low 
temperatures (T < 10^ K), and it becomes the dominant mech- 
anism for the damping of r-modes of cooler stars. It results 
from the momentum transport caused by the particle-particle 
scattering. In general several scattering processes can con- 
tribute to the total shear viscosity which can be approximately 
written as a sum of the partial shear viscosities of each individ- 
ual process. It has been widely thought that in a normal fluid 
star shear viscosity is completely dominated by the neutron 
scattering. Therefore, in the analysis of the r-mode instability 
in nonsuperfluid stars, it has been usual to take r] = T7„, using 
for rj„ the one calculated by Flowers and Itoth lf86l [stIi . and 
fitted by Cutler and Lindblom |88] in the simple form 
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FIG. 2: (Color online) Bulk (panel a) and shear (panel b) viscosities as a function of the symmetry energy slope parameter L for several 
densities and different models. Solid lines show the power-laws S, = A^lPi and r\ = AyjlP^ (see the text). The frequency of the mode and the 
temperature are taken 10^ s^' and 10^ K, respectively. In the inset it is shown the lepton fraction as a function of L for the same densities and 
models. 



where p\s and Tg mean that the density and temperature are 
given in units of lO'^ g cm^^ and 10^ K, respectively. Never- 
theless, recently, Shternin and Yakovlev [89.1 have shown that 
the main contribution to the shear viscosity at the temperatures 
relevant for the spin-down evolution of neutron stars comes 
from electron scattering when Landau damping is taken into 
account in the collision of charged particles mediated by the 
exchange of transverse plasmons. This electron contribution 
can be written as 1 18, 89] 

Tj, = 4 X 10-26(xp«fo) 14/97-5/3 g cm-1 s-i , (5) 

with Xp being the proton fraction, «i the baryon number den- 
sity in units of cm^^, and the temperature given in kelvins 
. Note that the temperature dependence of this contribution 
differs from the standard Fermi-liquid one Tje ^ T^^. In this 
work we take into account both contributions, Tj„ and Tjg. 

Fig.[T|shows the bulk (panel a) and shear (panel b) viscosi- 
ties for the BHF calculation of nonsuperfluid j3-stable npepL 
matter as a function of the density for T = 10^ K and a fre- 
quency of the mode o) = lO'* s^^ The contributions to the 
bulk viscosity from MURCA and DURCA processes involv- 
ing electrons and muons, as well as the contributions to the 
shear viscosity from neutron and electron scattering, are in- 
cluded. At low densities the only contribution to the bulk vis- 
cosity is due to the electron MURCA processes. The appear- 
ance of muons at ^ 0. 14 fm^^ switches on the muon MURCA 
processes with the consequent increase in the bulk viscosity. 
Electron and muon DURCA processes open in a jump-like 
manner at higher densities (^ 0.5 and ^ 0.59 fm^^, respec- 
tively for this model) when Xp > xourca- The main contri- 
bution to the shear viscosity, as it was said, comes from the 
electron scattering which is only exceed by Tj„ at densities 
rih < 0.17 — 0.18 fm^^. We note that, although rje dominates 
in general over )]„, the contribution of the neutron scattering 
can be larger than that of the electron one for temperatures 
r < 10^ K, and, as it is seen in the figure, for subsaturation 
densities \M\- The results for APR, AFDMC, the Skyrme 



forces, and the RMF models are qualitatively similar to the 
ones obtained within the BHF approach and, therefore, we do 
not show them here for simplicity. 



In Fig. |2] we show the dependence of the bulk (panel a) 
and shear (panel b) viscosities on the symmetry energy slope 
parameter L for densities 0.08,0.16 and 0.32 fm^^, and the 
different models considered. As in Fig. [T] the frequency of 
the mode and the temperature are taken lO'* s"' and lO'' K, 
respectively. The figure shows that both viscosities increase 
with L for all densities, except for the lowest one for which 
the shear viscosity decreases. This is just a consequence of the 
dependence of both viscosities on the lepton fraction (see Eq. 
©, Eq. (35)of Ref. [84], andEqs. (18-21) of Ref. [85]) which 
increases with L above saturation density, and decreases be- 
low it, as it can be seen in the inset of the figure. The depen- 
dence of £, and Tj with L can be described by simple power- 
laws of the type ^ — A^L^^ and rj — Aj^lP^ at each density, 
shown by solid lines in the figure. For completeness, we plot 
in Fig. [3]the density dependence of the coefficients A^ and A-q 
(panel a) and the exponents and Brj (panel b). The contri- 
butions to t, from MURCA and DURCA processes are shown 
separately. Note that A^^^^^^ and B^^^^^^ are only defined 
for densities larger than the DURCA threshold. It is observed 
that A^Pj^^^^ increases until it reacheas a plateau at '-^ 1 fm^^, 
while Bf„,,„^. always decreases. On the other hand, Ap,,,,„^, 

^DURCA -' ^MURCA 

^^^MURCA^ decreases (increases) initially, then reaches a mini- 
mum (maximum) around ~ 0.45 fm^^, and finally increases 
(decreases). A similar behavior is observed for A,j and Bjj. 
Howerver, note that B,j is negative below saturation density, 
contrary to B^^^^^^ and B^^^^^^ which are always positive. 
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FIG. 3: (Color online) Density dependence of the coefficients and Ajj (panel a) and the exponents and Sjj (panel b). The circles and 
squares show, respectively, the results of the modified (^murca) and direct (^durca) Urea contributions to whereas the triangles display 
those of 7] . As in the previous figures the frequency of the mode and the temperature are taken CO = 10"^ s"' and T = lO'^ K. 



IV. DISSIPATIVE TIME SCALES OF R-MODES 



The dissipative time scale of an r-mode is given by |@] 



2E \dt J 



(6) 



where the index / refers to the various dissipation mechanisms, 
in our case bulk viscosity, shear viscosity and gravitational 
wave emission, E is the energy of the mode, and {dE/dt)i is 
the rate of dissipation associated with each mechanism. The 
energy E can be expressed as an integral of the fluid perturba- 
tions 1i6i,lZi] 



p5v-8v* + 



5p 



50 5p' 



d^r 



(7) 



with p being the mass density profile of the star, and 5v, 5p, 
5<I> and 5p the perturbations of the velocity, pressure, gravita- 
tional potential and density due to the oscillation of the mode. 
For the case of r-modes in the small angular velocity limit E 
can be reduced to a simple one-dimensional integral 16|], 



2 Jo 



(8) 



where a is the dimensionless amplitude of the mode, and R 
and n are the radius and the angular velocity of the star, re- 
spectively. Here we focus only on r-modes with angular quan- 
tum number / = 2 and azimuthal projection m = 2 since, as 
shown e.g., in Refs. 10, r-modes with I = m = 2 are the 
dominant ones. Higher multipoles lead to weaker instabilities, 
and are not considered in this work. 

The dissipation rate due to the bulk viscosity is given by 



dE\ 
'dfj 



^|V.5v|Vr 



(9) 



In general, the quantity | V • 5v'p is a complicated function of 
the radial and angular coordinates. However, for slow rotating 



stars, the bulk viscosity ^ depends to lowest order only on the 
radial coordinate. Therefore, it is usual (see Refs. |IB.l3.H>[i2l) 
to define the angle-average (|V ■ dv\'^) which allows also to 
reduce [dE/dt)^ to a one-dimensional integral. 



dE\ 
dt ) I 



An I ^{\V-8v\^)r^dr 



(10) 



The quantity (| V • 5v'p) can be determined numerically ||7tl- In 
this work, however, we will use the analytic expression given 
by Lindblom and Owen in Refs. ll9[ |T2ll . 



{\V.5v\^)^^(!^y 1+0.86^ 
^' ' ' 690 \rJ \rJ 



\nGpJ 



(11) 

In this expression, p = M/{4kR^/3) is the average density 
of the nonrotating star, and G is the gravitational constant. 
Finally, using Eqs. (|8]l, (fTOl i and ( fTTT l we get 



4^ ( \\21- 
690 V ^Gp J 







pr^'+^dr 



-1 



l+0.86(^'' 



r^dr . 



(12) 



The dissipation rate due to the shear viscosity is given by 



dE 
dt 



2 I n5a"'5a:,d\ , 



where the SOah is the shear defined as tSl 



1 / 2 



(13) 



(14) 



Working out the angular integrals (see Refs. ||l,[T3l[T3l) and 
using Eq. ([8]l yields 



— = (/-l)(2/ + l) 



pr^'+^dr 



/ r]r^'dr. (15) 
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FIG. 4: (Color online) Dissipative time scales as a function of the temperature for a IAMq (panel a) and 2M0 (panel b) neutron star rotating 
at 10% of its Kepler frequency. Results are shown for three Skyrme forces. Damping time scale due to bulk (shear) viscosity is shown by the 
solid (dashed) lines. The time scale associated with the groth of the mode due to the emission of gravitational waves Tqw is shown by the 
horizontal dotted lines. The frequency of the mode is taken CO = 10"* s^' . In parenthesis it is given the value of the slope parameter L of each 
model. 



Finally, the time scale of the growth of an r-mode due to 
the emission of gravitational waves is given by |01 



1 



32;rGn 



21+2 



(/-I) 



r.2/+3 



[(2Z + 1)!!]2 + l 



1 + 2 



2/+1 



(16) 

We plot in Fig. |4] the time scales , T,] and Tqw as a 
function of the temperature for a 1 AMq (panel a) and 2Mpj 
(panel b) neutron star rotating at 10% of its Kepler frequency 



{Q.Kepier ~ 1 SOO {M /Mq) (lOkm/ R)^ S-' IMS])- As an 
example, we show results for the Skyrme forces SLy 10, SkMP 
and SkT4 which give values of L = 39.2,69.7 and 93.4 MeV, 
respectively. Other Skyrme forces, RMF models, and the 
BHF, APR and AFDMC calculations give similar qualitatively 
results, and are not shown here for simplicity. As in the pre- 
vious figures the frequency of the r-mode is taken 10^ s^'. 
Note first that Tgw is larger for the models which give a larger 
value of L. This is because a larger value of L implies a stiffer 
EoS and, therefore, a less compact neutron star (i.e., a more 
extended and less dense object). Consequently (see Eq. (fTSIl), 
l/tQw is smaller and Zqw is larger According to Eqs. (fT2] i 
and ( fTsT i and Ttj decrease when increasing £, and Tj, respec- 
tively. However, we have just seen that ^ and rj increase with 
L so, contrary to Tgw^ the models with larger L predict smaller 
values of and Tj; . Finally, we note that the three time scales 
decrease when increasing the mass of the object. In fact, for a 
given EoS, the more massive is the star, the denser it is. Then, 
it is clear from Eq. (fTSI l that Tgw decreases. Morever, £, and rj 
increase also with the mass of the star because of their increase 
with density (see Fig. [T]|. Assuming constant profiles for the 
density, ^ and Tj, one can see from Eqs. (fT2l) and (fTSl l that 
and T,j behave as ^{p/^ )R'^ and T,j '^{p/7] )R^. Since the 
increase of ^ and r\ with the mass of the star is much faster 
than that of p this explains the decrease of and i-q ■ 



V. R-MODE INSTABILITY REGION 

The time dependence of an r-mode oscillation is given by 
^im-t/x ^ where CO is the frequency of the mode, and T is an 
overall time scale of the mode which describes both its expo- 
nential growth, driven by the CFS mechanism ll20ll2lll . and its 
decay due to viscous damping ^ . It can be written as 



1 



1 



1 



+ ■ 



1 



T(i2,r) TGw(ii) T^(n,r) T^(r 



(17) 



If Tgw is shorter than both and the mode will expo- 
nentially grow , whereas in the opposite case it will be quickly 
damped away. Therefore, it is clear that the r-mode will be 
stable only when 1/t is positive. For each star at a given tem- 
perature T we can define a critical angular velocity £2^ as the 
smallest root of the equation \/t{Q.c,T) ~ 0. This equation 
defines the boundary of the so-called r-mode instability re- 
gion. A star will be stable against the r-mode instability if its 
angular velocity is smaller than its corresponding £2^ . On the 
contrary, a star with £2 > £2^ will develop an instability that 
will cause a rapid loss of angular momentum through gravita- 
tional radiation until its angular velocity falls below the criti- 
cal value. 

In Fig. |5] we present the r-mode instability region for a 
1.4M0 (panel a) and 2Mq (panel b) neutron star obtained 
for some Skyrme forces (solid lines), RMF models (dashed 
lines), and the BHF (dotted line), APR (dotted-dashed line) 
and AFDMC (double-do tted-dashed line) calculations. We 
note that the BHF and APR results are not shown for the 2Mq 
neutron star because the maximum mass predicted by these 
models is ^ LSM;; and ^ \.92Mp„ respectively. We note 
also that the value 1 .92Mq is sligh tly lower than the 2.2Mq of 
the original APR calculation 11461] . the reason being, as it was 
already mentioned, that in the Heiselberg and Hjorth-Jensen 
parametrization |58] the EoS is softened at higher densities 
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Temperature T [K] Temperature T [K] 



FIG. 5: (Color online) r-mode instability region for a 1.4M0 (panel a) and IM-;., (panel b) neutron star obtained for some Skyrme forces 
(solid lines), RMF models (dashed lines), and the BHF (dotted line), APR (dotted-dashed line), and AFDMC (double-dotted-dashed line) 
calculations. The frequency of the mode is taken CO = 10"* s^^ . In parenthesis it is given the value of the slope parameter L of each model. 



in order to obey causality. The angular velocity is given in 
units of Hfcepieri ^nd the value of the slope parameter L of 
each model it is shown in parenthesis. The r-mode instability 
region is larger for the more massive star. This can be un- 
derstood from our previous discussion of Fig. |4] Note there 
that, with i2 fixed, the crossing point between Tgiv and 
(Tj)) moves to higher (lower) temperatures when going from a 
1 AMqj to a 2M(;) star, therefore, making the instability region 
wider. Dissipation due to shear viscosity kills the mode at low 
temperatures, while the bulk viscosity does it at high ones. In 
fact, shear viscosity supresses completely the r-mode insta- 
bility for temperatures below 10^ K. Similarly, bulk viscosity 
prevents the mode from growing in a star that is hotter than 
a few times lO'" K. For temperatures between these two the 
growth time due to gravitational radiation is short enough to 
overcome the viscous damping, and drive the r-mode unsta- 
ble. Note that the instability region is smaller for the models 
which give larger values of L. The reason is simply the fact 
that both bulk and shear viscosities, as we already discussed, 
increase with L and, therefore, the damping of the mode is 
more efficient for the models with larger values of L. 

Finally, we plot in Fig. |6] the dependence on L of the crit- 
ical angular velocity for a fixed temperature, taken equal to 
the estimated core temperature of the pulsar in the LMXB 4U 
1608-52 (hereafter called simply 4U 1608-52), T -4.55 x 10^ 
K [94]. Results of the different models are shown for two 
possible values of the mass of this object, 1.4Mq (panel a) 
and 2Mq (panel b). The horizontal lines show the observa- 
tional spin frequency of 4U 1608-52 (620 Hz 1,95.1) in units of 
^Kepler, assuming that its radius is: i) 10, ii) 1 1 .5, iii) 12 or iv) 
1 3 km. Most of the rapidly rotating neutron stars in LMXB are 
observed to rotate at spin rates well below Sleepier- Although 
some of them can reach spin frequencies larger than £2^, it 
is expected that they spend a very short time inside the in- 
stability region since they would rapidly spin down due to the 
emission of gravitational waves lll7[l . Therefore, most of these 
objects should likely be outside the instability region li96ll97ll . 



It is clear from the picture, then, that if the radius of 4U 1608- 
52 is smaller than '--^ 1 1 .5 km, this object is always out of the 
instability region for any model (£2c is larger than its spin fre- 
quency), and we cannot conclude anything about the value 
of L. On the other side, if its radius is larger than ^ 12(13) 
km and its mass 1.4M0(2M0), 4U 1608-52 is always inside 
the instability region {£lc is smaller than its spin frequency), 
and we can neither draw any conclusion on L from this pul- 
sar. Only if its radius is in the range 1 1 .5 — 12(1 1 .5 — 13) km 
and its mass 1.4M0(2M0) we can say that observational data 
seem to favor values of L larger than ^ 50 MeV, if 4U 1608- 
52 is assumed to be outside the instability region. This is in 
contrast with the recent work of Wen, Newton and Li llisll 
where they show, as we said, that smaller values of L seem 
to be more compatible with observation. We should mention, 
however, that these authors assume that the main dissipation 
mechanism of the r-mode is due to the viscous boundary layer 
at the crust-core interface where densities are smaller than po. 
Therefore, their calculation of the shear viscosity is done in 
a region of densities for which, as we saw, T] decreases with 
L. Consequently, they obtain that r-mode instability region is 
smaller for smaller values of L. Nevertheless, in our work, 
we need to calculate the bulk and shear viscosities in a range 
of densities covering the whole density profile of the star to 
determine their corresponding damping time scales (see Eqs. 
(fT2l i and (fTSTi). Both viscosities increase with L in this range 
of densities and, therefore, we reach a conclusion opposite to 
that of the authors of Ref. llisll . 



VI. SUMMARY 

In this work, we have studied the role of the sym- 
metry energy slope parameter L on the r-mode instabil- 
ity of neutron stars. To such end, we have used differ- 
ent models for the nuclear EoS that include the microscopic 
Brueckner-Hartree-Fock approach, the variational Akmal- 
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FIG. 6: (Color online) Critical angular velocity as a function of the symmetry energy slope parameter L for a XAMq (panel a) and 2M0 (panel 
b) neutron star at the estimated core temperature of 4U 1608-52, T ~ 4.55 x 10^ K |94], and different models. The frequency of the mode is 
taken ft) = 10'* s^'. Solid lines show the result of a quadratic fit. The horizontal dashed-lines show the observational spin frequency of 4U 
1608-52 in units of O-Kepley assuming that the radius of this object is: i) 10, ii) 1 1 .5, iii) 12 or iv) 13 km. 



Pandharipande-Ravenhall EoS, a parametrization of recent 
Auxiliary Field Diffusion Monte Carlo calculations, and sev- 
eral phenomenological Skyrme forces and relativistic mean 
field models. We have found that the r-mode instability re- 
gion is smaller for those models which give larger values of 
L. We have shown that this is due to the fact that both bulk 
and shear viscosities increase with L and, therefore, make the 
damping of the mode more efficient for the models with larger 
L. We have shown also that the dependence of both viscosi- 
ties on L can be described at each density by simple power- 
laws of the type E, ~ A^L^^ and Tj Aj^lP^ . Finally, we have 
tried to constrain the value of L using the measured spin fre- 
quency and the estimated core temperature of the pulsar in the 
low-mass X-ray binary 4U 1608-52. We have concluded that 
observational data seem to favor values of L larger than ^ 50 
MeV if this object is assumed to be outside the instability re- 
gion, its radius is in the range 11.5 — 12(11.5 — 13) km, and 
its mass \AMc.j{2Mqj). Outside this range it is not possible to 
draw any conclusion on L from this pulsar. These results are 



in contrast with the recent work of Wen, Newton and Li ll45ll . 
where these authors show that observation seems to be more 
compatible with smaller values of L. Finally, we note that 
the inclusion of o ther sourc es of dissip ation, such as e.g., hy- 
peron fllMJTol or quark ifToelfrTell bulk viscosities, is not 
expected to change the qualitative conclusions of this work. 
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